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High resolution, inelastic neutron scattering measurements on SrCu2 (603)2 reveal the dispersion 
of the three single triplet excitations continuously across the (H,0) direction within its tetragonal 
basal plane. These measurements also show distinct Q dependencies for the single and multiple 
triplet excitations, and that these excitations are largely dispersionless perpendicular to this plane. 
The temperature dependence of the intensities of these excitations is well described as the comple- 
ment of the dc-susceptibility of SrCu2 (603)2 . 
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Quantum magnets which display collective singlet 
ground states have been of much recent interest^. Sev- 
eral of these, such as CuGeOgS, MEM-(TCNQ)2^, and 
NaV205^ result from spin-Peierls phenomena in low di- 
mensions where the lattice combines with s=l/2 spin de- 
grees of freedom to break translational symmetry below 
some characteristic phase transition temperature, and a 
non-magnetic ground state with a characteristic energy 
gap is observed. A related state is observed in the s=l 
antiferromagnetic chain compounds, such as CsNiCla^ 
and NENP— , where a non-magnetic ground state with 
an energy gap, the Haldane gap, forms in the absence of 
translational symmetry breaking. 

SrCu2 (603)2 has been proposed2*^ as a realization of 
the Shastry-Sutherland model^ of interacting dimers in 
two dimensions. This material crystallizes^" into the 
tetragonal space group I42m with lattice parameters 
a=8.995 A, c=6.649 A. Magnetically, it can be thought 
of in terms of well isolated basal planes populated by an- 
tiferromagnetically coupled s=l/2 moments on the Cu^+ 
sites. These are arranged in dimers at right angles to 
each other and forming a square lattice. The micro- 
scopic Hamiltonian appropriate to this system is based 
on: 
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It is known that both interactions, J and J', are anti- 
ferromagnetic and similar in strength, such that the sys- 
tem is not far removed from the critical value of x=J'/J 
known to be appropriate to the quantum phase tran- 
sition which separates a 4 sublattice Neel state from a 
collective singlet ground state. 

The estimated values of J and J' have evolved over 
time as both theory and experiment have improved. 
Early estimates of J and x were J=8.6 meV and x=0.68i, 
very close to the critical value Xc=0.69^^, where the sin- 
gle triplet excitation goes soft. More recently Miyahara 
and Uedcki^ found J=7.3 meV and x=0.635 from fits to 



the magnetic susceptibility. Somewhat smaller values 
were obtained by Knetter et aU^ who compared theoret- 
ical and experimental ratios of the energy of the lowest 
S=l two triplet bound state to the single triplet gap, 
to obtain J=6.16 meV and x=0.603. Note that within 
this theorjii^, the collective singlet ground state becomes 
unstable when the lowest energy two triplet bound .state 
goes soft. 

Earlier, relatively low resolution inelastic neutron 
scattering measurements^'^'^^'^^have directly identified 
three bands of excitations corresponding to single (n=l) 
triplet excitations, as well as to two (n=2), and to three 
(n=3) triplet excitations. These measurements directly 
show the appearance of the energy gap in the spectrum of 
excitations with decreasing temperature, as well as the 
dispersion of these excitations in an applied magnetic 
field. 

The more recent of these studieai^tiSi has been able 
to investigate subleading terms in the spin Hamiltonian. 
Terms such as the Dzyaloshinski-Moriya (DM) interac- 
tion, which is allowed by symmetry between spins on 
neighboring dimers, and anisotropic exchange, can be re- 
sponsible for both dispersion of these excitations, and the 
removal of the three-fold degeneracy which would oth- 
erwise characterize the n=l triplet excitation spectrum 
in SrCu2 (603)2 ■ The presence of such small terms in 
the spin Hamiltonian has also been investigated through 
high field specific heat measurementsii, performed on 
samples from the same crystal growth as that used in 
the present study, and through ESR measurementsiSiSS. 

In this letter, we report new high resolution inelastic 
neutron scattering measurements, which probe both the 
energy and Q dependencies of these previously identified 
bands of excitations with new precision. Measurements 
were performed on two different high resolution cold neu- 
tron instruments, allowing both sufficiently high energy 
resolution to resolve the three n=l triplet excitations 
in SrCu2 (603)2 , and sufhciently high Q resolution to 
discern different Q-dependencies of the n-triplet excita- 
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tions, where n=l,2,3. 

The present smgle crystal of SrCu2(^^B03)2 was 
grown from a self-flux by floating zone image furnace 
techniques in an O2 atomosphere using a four mirror 
furnace. It is cylindrical in shape, with approximate 
dimensions 0.6 cm in diameter by 10 cm long, and an 
a- axis was roughly aligned along the cylindrical axis. 
Small pieces of the crystal were used for bulk charac- 
terization, and the characteristic fall off of the dc sus- 
ceptibility near 10 K, as well as well defined steps in the 
magnetization versus applied magnetic field at high field, 
were observedii. Neutron diffraction measurements, en- 
abled by the use of ^-^B isotope, revealed a high quality 
single crystal throughout the volume of the sample with 
a mosaic spread of less than 0.2 degrees. 

The crystal was mounted in a pumped ^He cryo- 
stat with the long cylindrical axis vertical, placing the 
(H,0,L) plane of the crystal coincident with the hori- 
zontal plane, and hence the scattering plane. Neutron 
scattering measurements were performed using the Disk 
Chopper Spectrometer (DCS) and the SPINS triple axis 
spectrometer, both located on cold neutron guides at the 
NIST Center for Neutron Research. 

The DCS uses choppers to create pulses of monochro- 
matic neutrons whose energy transfers on scattering are 
determined from their arrival times in the instrument's 
913 detectors located at scattering angles from -30 to 
140 degrees. Using 5.1 meV incident neutrons and the 
medium resolution optioni^, the energy resolution is 
0.09 meV. 

The SPINS triple axis spectrometer was operated us- 
ing seven pyrolitic graphite analyser blades accepting 7 
degrees in scattering angle, and neutrons of 5 meV, fixed 
scattered energy. A cooled Be filter was placed in the 
scattered beam to remove contamination of higher order 
neutrons, and the resulting energy resolution was ^ 0.5 
meV. 

A program of constant-Q measurements across the 
(11,0, L) plane of SrCu2(B03)2 at 1.4 K was carried out 
covering energies from 1 to 8 meV. A compendium of 
such scans with intervals of AH=0.2 is shown in the color 
contour map in Fig. la, which displays data within the 
tiuj, H plane at L=0. Repesentative scans making up 
this map are shown in Fig. lb. One clearly identifies 
both the n=l triplet excitation near 7ia;=3.0 meV and 
the n=2 triplet excitation, whose spectral weight is max- 
imum near 4.9 meV. One also sees a continuous 
component to the n=:2 triplet excitation, which extends 
to the highest energies collected in this set of measur- 
ments, 8 meV. 

These measurements are qualitatively similar to those 
first reported by Kageyama et al>i^, within the (II,K,0) 
plane of SrCu2(B03)2 and with lower energy resolution; 
however there are important differences. For example, 
as can be seen in Fig. 1, we observe no substantial dis- 
persion of the maximum of the spectral weight of the 
n=2 triplet excitation in contrast to a bandwidth of 1.5 
meV reported in these earlier measurements^"'. Also, as 
we discuss next, the Q-dependence of these excitations 
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FIG. 1: Top panel: A map of the measured dynamic struc- 
ture factor for SrCu2 (603)2 at T=1.4 K along the (H,0,0) 
direction is shown. The intensity scale is logarithimic. The 
map is made up of constant-Q scans of the form shown in 
the bottom panel, where Q=(1.5,0,0). 

is different from that reported by Kageyama et aU^. 

Figure 2 shows the results of constant energy scans 
performed at huj=3, 4.8 and 9 meV, corresponding to the 
n-triplet excitations with n=l,2,3 respectively. These 
measurements clearly show the n=l excitation at 3 meV 
to peak up at half integer values of H, that is at 11=1.5 
and 2.5, in maps of this scattering within the (H,0,L) 
plane, and at integer values of 11=2 for the n=2 excita- 
tion at 4.8 meV and for the n=3 excitation at 9 meV. 
These results indicate distinct form factors for the n- 
triplet excitations, with the n=:l triplet being different 
from the multi-triplet excitations. Knetter and Uhrig^i 
have recently calculated the n=2 triplet contribution to 
the dynamic structure factor within the (II,K,0) plane 
of SrCu2(B03)2 by perturbative techniques. They show 
that it is peaked at (2,0,0), consistent with the present 
set of inelastic scattering measurements. 

It is also clear from Fig. 2 that the Q-dependence of 
all the excitations show little L dependence, consistent 
with well isolated two dimensional basal planes. Given 
that the scattering displays very little L-dependence, the 
DCS measurements can be integrated along L, resulting 
in a good quality determination of the dispersion of the 
n=l triplet excitations in the (H,0) direction within the 
tetragonal basal plane. This is what is shown in Fig. 
3, where the top panel shows a color contour map of 
the inelastic scattering acquired using 5.1 meV incident 
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FIG. 2: Constant energy scans at 3 (top panel), 4.8 meV 
(middle panel) and 9 meV (bottom panel) probing the 
Q-dependence of the n=l,2, and 3 triplet excitations in 
SrCu2(B03)2 within the (H,0,L) plane at T=1.4 K. 



neutrons. The bottom three panels show cuts through 
this map, which approximate constant-Q scans at (-1,0), 
(-1.5,0) and (-2,0), from top to bottom, respectively. 

These inelastic measurements clearly resolve three 
branches of excitations, corresponding to S^=0,±1 tran- 
sitions from the ground state, the lower two of which 
are nearly degenerate at (-1.5,0). These measurements 
have the advantage that they track the dispersion of the 
single singlet-triplet excitations continuously in Q along 
(H,0). This color contour map is also consistent with 
the mapping of the Q-dependcnce of the integral of the 
singlet-triplet scattering, shown in Fig. la, which shows 
this scattering to peak in intensity near (-1.5,0). 

Although the high resolution n— 1 triplet dispersion 
curves shown in Fig. 3 have not been previously re- 
ported, a substantial amount of information is known 
about these spectra from ESR measurements. In partic- 
ular Nojiri et al.^" used ESR to measure the q=0 n=l 
triplet excitations finding states at 2.81 meV (679 ± 2 
GHz) and 3.16 meV (764 ± 2 GHz), in excellent agree- 
ment with the lowest and highest n=l triplet excitations 
at Q=(-1,0) and (-2,0) shown in Fig. 3. The splitting be- 
tween these two excitations has been shown to be 4Di^, 
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FIG. 3: Top panel; Colour contour map of the dynamic struc- 
ture factor for the n=l triplet excitations along the (H,0) di- 
rection within the basal plane of SrCu2 (603)2 . These mea- 
surements, taken with the DCS spectrometer, integrate along 
L, as described in the text. Bottom panel: Cuts through the 
map of the top panel are shown, which approximate constant- 
bf Q scans and clearly resolve the three branches to the n=l 
triplet excitation. The solid lines are Gaussian fits to the 
data. 



where D is the inter-dimer DM interaction. In fact the 
dispersion relation shown in Fig. 3 bears similarities to 
the calculation shown in Fig. 2c of Cepas et al>i^, except 
for the small gap, ^ 0.2 meV, we observe between the 
S^=0 transition and the upper of the S^=±l branches 
at Q=(-1.5,0). This gap must result from some other 
anisotropic interaction, such as the recently proposediS 
perpendicular intra-dimer DM interaction. 

We note that the band width of the middle, S^=0 
mode in Fig. 3 is extremely small, roughly 0.05 meV. 
This mode is the least affected by anisotropic interac- 
tions, and hence it is the mode most directly compara- 
ble to calculations of the n=l triplet excitation based on 
Eq(l). This bandwidth, which scales like x^, is similar 
to the bandwidth found by Weihong et al,^^ for x=0.6, 
which is about 4% of the n=l triplet gap energy. 

We have also measured the temperature dependence 
of both the n=l and n=2 triplet excitations with the 
SPINS triple axis spectrometer, and that of the n=l 
triplet excitations with DCS. The SPINS measurements 
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FIG. 4: The temperature dependence of inelastic intensity at 
(1.5,0,0) and fiu) =3 meV, as well as at (2,0,0) and ?icij=4.8 
meV are shown. The inset shows the raw intensity data, 
while the main figure shows the normalized intensity assum- 
ing zero at 25 K. This is compared to the complement of the 
measured dc susceptibility (l-x)i which clearly describes well 
the observed temperature dependence. 

at Q=(1.5,0,0) and (2,0,0) and energy transfers of Ticu 
—3 meV and 4.9 meV, sample the n=l and n—2 triplet 
excitatons, respectively, and are shown in Fig. 4. The 
DCS measurements integrate the inelastic scattering in 
data sets of the form shown in the top panel of Fig. 3. 



The DCS measurements integrate between 2.7 and 3.3 
meV and across all wavevectors from H=-2.25 to H=- 
0.75 along (H,0) within the basal plane. 

Figure 4 shows that the temperature dependence 
of the n=l and n=2 triplet excitations are identi- 
cal, which may not have been concluded from earlier 
measurements^^. We do however confirm the very rapid 
dropoff of this inelastic intensity with increasing tem- 
perature. Such a strong temperature dependence is not 
characteristic of a system undergoing a phase transition, 
and indeed there is no evidence of a broken symmetry 
associated with the appearance of the collective singlet 
ground state. Rather the temperature dependence of 
the inelastic scattering can be very well described as the 
complement of the dc susceptibility. That is that the 
temperature dependence to the inelastic scattering fol- 
lows that of the dc susceptibility. This has been mea- 
sured and modelled using finite temperature Lanczos 
methodii. The complement of %, referred to as 1-x in 
Fig. 4, is given by x(T=20 K)-x(T). This quantity is 
scaled to compare to the temperature dependence of the 
inelastic scattering, and one can see that this provides 
an excellent description of the temperature dependence 
of the inelastic scattering. 
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